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Abstract 

Frequency-domain nonlinear wave mixing processes may be described either using response func- 
tions whereby the signal is generated after all interactions with the incoming fields, or in terms 
of scattering amplitudes where all fields are treated symetrically with no specific time ordering. 
Closed Green's function expressions derived for the two types of signals have different analytical 
properties. The recent controversy regarding the sign of radiative damping in the linear (Kramers 
Heisenberg) formula is put in a broader context. 
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I. INTRODUCTION 



A lively debate is currently going on with regard to the correct sign that should be in- 
troduced in the Kramers-Heisenberg expression for the (linear) optical polarizability. Both 
"opposite sign" and a "constant sign" formulas have been derived by various authors. Per- 
turbative QED calculations have been carried out to include radiative damping [1-8]. In 
an insightful recent article Bialyncki-Birula and Sowinski [9] had pointed out that this issue 
is fundamentally connected with the linear response vs. the scattering points of view for 
Rayleigh scattering. In this letter we extend this argument to nonlinear wave mixing of ar- 
bitrary order. Closed formal expressions are derived that reveal the analytical properties of 
both retarded (response) and non retarded (scattering) signals. The linear response results 
are recovered to lowest order. Our microscopic deviation could serve as the starting point 
for a full QED calculation. However, this is not required in order to pinpoint the analytical 
properties which should hold for other types of nonradiative damping as well. 

The following derivation applies for processes of arbitrary order. However, for clarity we 
focus on four wave mixing. Consider a system (atom, molecule) interacting with four modes 

of the radiation field. The states of the system will be denoted a,b,c The j'th mode 

of the radiation field has a frequency uij and initial occupation number rij. To simplify the 
notation we hereafter consider sum frequency generation whereby = uji+uj2+uj 3 , however 
the results can be easily extended to any combination of frequencies u; 4 = ±a>i ± u> 2 + ±u; 3 . 
For the process of interest the initial state of the system + field is | i >—\ a, nxn^nzn^ >, 
and the final state is | / >=| a, ri\ — l,n 2 — l,n 3 — l,n 4 + 1 >, with energies E { and Ef 
respectively. 

The coupling of an atom located at point r with the radiation field is 



We shall divide the electric field into its positive and negative frequency components 



H mt = -V.E(r,t). 



(1) 




(2) 




(3) 



2 



cij is the photon annihilation operator of mode j, V is the dipole operator and Q is the 
quantization volume. 

We first assume that the system is prepared in a nonequilibrium steady state with the 
four fields and consider the entire process as a single four-photon scattering event. This 
is a non causal process whereby all four modes are treated along the same footing. This 
implies that the interaction with model 4 need not be the last, and all time orderings are 
allowed and should be summed over. The process is described by the S matrix element 

S fi = A f sr fi (E i )8(E i -E f ), (4) 

where T/j is the matrix element of the T matrix T=V+V G(E) V, and 

is the retarded Green's function. The scattering amplitude for this process is given by 

4? = ^E ( 6 ) 

pi 



< a | VG(E a + u! + uj 2 + uj 3 )VG(E a + wi + u 2 )VG(E a + ui)V | a > %i + u 2 + u 3 - co 4 ) 
where 

Afi = (~^~) 2 V n l n 2^3(n4 + 1) V / ^1^2^3^4- (7) 

^ denotes the sum over all 4! permutations of u>i, cu 2 , <^3, and -oj±. The sign convention 
for Uj in eq.6 is as follows: an absorbed photon gives +Uj whereas an emitted photon gives 
—ujj. In the process considered here, wi,o>2 and uj 3 are absorbed and o> 4 is emitted. Other 
processes can be calculated by simply changing the signs of a;j,as warranted in each case. 
H is the Hamiltonian for the atom + all modes of the radiation field excluding the 4 modes 
of interest, since those were treated perturbatively. 

We next turn to the standard semiclassical description of four wave mixing whereby the 
system first interacts with modes 1, 2 and 3 to generate the third-order polarization which 
then serves as a source for mode 4 [10]. This is a causal response where mode 4 is special 
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since its interaction with the system must be the very last. We shall calculate the n'th 
order polarization as the expectation value of V [11] 

n 

pin) _ < \p( m ) | V | ^ (n ~ m) >, (8) 

m=0 

where ^( m ) is the perturbed wavefunction to m'th order in V. The m'th term represents 
(n — m)'th order for the ket and m'th order for the bra. Overall there are (n + 1) terms. 
For n = 3 we get four terms corresponding to m — 0, 1, 2, 3 respectively in eq.(8). We then 
get = Afi xfi where 



P3 

< a | VG(E a + ui + uj 2 + uj 3 )VG{E a + uj x + u 2 )VG{E a + wi)V | a > 

< a | VG\E a + ^+002- uj 4 )VG{E a + wi + cj 2 )\/G(£ a + u;i)V | a > 
+ < a | FG f (wi - u 4 + u 2 )VG\E a + uj 1 - u 4 )VG(E a + u x )V \ a > 

)VG\E a - u 4 + u{)VG\E a - uo A )V I a> 



is the susceptibility and 



is the advanced Green's function. J2 P3 denotes the sum over all 3! permutations of the 
incoming modes 001,002 and 003. In this expression the sign conversion of 00 j is the same 
as in eq. (6) (+sign for absorbed photons, - for emitted). Each coupling with the ket is 
accompanied by a retarded Green's function G, whereas an advanced Green's function G^ 
is accompanied by bra interactions. Eq. (8) will thus yield (n — m) G and m G^ factors. 
An analogous expression was obtained recently using superoperators in Liouville Space. 
[12] The present Hilbert Space form is more suitable for comparison with the scattering 
amplitudes. 

The extension of Eqs.(6) and (9) to arbitrary order is straightforward. 1 5^ n+1 - ) will have a 
single basic term with (n + 1) V factors, and n retarded Green's functions G with arguments 
E a +coi, ...E a +uj\+uj2 + ■■■+to n ,. In addition, it contains a sum over the (n+1)! permutations 
of coi,..., oo n +i- Each tOj may be changed to — coj to describe different processes, nas 
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(n + 1) basic terms each containing (n + 1) V factors, and a (G + ) m (G) n ~ m factor with 
m — 0, ...n (See Eq.8). The retarded Green's functions will only depend on the incoming 
frequencies Ui r ..,u n . All advanced Green's functions also depend on the signal frequency 
—u n+ i. Each basic term yields n! terms upon the permutation over u>i.. uj n . Altogether 
both S^ n+1 ^ and x^ contain (n + 1)! terms. For S^ n+1 ^ these come from the [n + 1)! 
permutations of the (n + 1) frequencies over a single term, x^ has (n + 1) basic terms, 
each containing n! permutations of the n "incoming" frequencies. Unlike x^ n \ 5 ,( - ri+1 - ) is 
symmetric with respect to all (n + 1) modes; the interaction with the signal field co n +i need 
not be chronologically the last. gO+i) only contains retarded Green's functions and all 
propagations are forward in time. The (n+ 1)! permutations take care of the possible time 
orderings of interactions with the various modes, x^ ? m contrast, depends on both retarded 
and advanced Green's functions which correspond to forward and backward propagations 
respectively along the Keldysh Schwinger loop [13-15]. 

Both expressions (6) or (9), can serve as a starting point for a full QED perturbative 
calculation, where the coupling with all modes of the electromagnetic field is included in 
the Hamiltonian H. This will result in damping terms. However, the analytical properties 
of the two signals are completely determined by eqs (6) or (9), and will be invariant to the 
level of approximation used for the radiative damping. When all frequencies are tuned off 
resonance, we can neglect the ±ie terms in the Green's functions, setting G — G\ The 
causal and non causal expressions are then identical. 

The ongoing damping sign controversy was restricted to the linear response of a two level 
system with ground state a and an excited state b [1-9]. This may be immediately resolved 
by our general formulation. For Rayleigh (elastic) scattering we have one absorbed and one 
emitted photon with the same frequency. We thus set u>i = uj and uj 2 = —UJ- The second 
order analogue of eq.(6) reads 



This is obtained from a single term + 2 permutations of uj and —to. Writing the matrix 
elements explicitly for the system Hamiltonian (neglecting coupling with other radiation 
modes), this gives 



,(2) 
'fi 



< a | VG{E a + oo)V \a> + <a \ VG(E a - uj)V \ a > 



(11) 




1 1 



(12) 



V ah 



E a + uj — E b + is E a — uj — E b + is 
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This is known as the "constant sign" prescription (both terms have a +is factor) [1-9]. 
The linear response expression similar to eq.(9) gives on the other hand 



This comes from two basic terms with no permutation. (there is only one incoming 
field which is absorbed and has frequency a;). Taking the matrix elements we recover the 
"opposite sign" prescription 



The origin of the different damping signs of and had been clearly pointed out 
by Bialynicki-Birulya and Sowinski [9] who had further carried out a fourth order QED 
calculation of dampling using both expressions. The present results extend these arguments 
to nonlinear processes of arbitrary order. It should be emphasized that the fundamental 
principle of causality must always hold and is never in doubt. The terms "causal" and 
"noncausal" in this letter refer to different observables. Traditional semiclassical formulation 
of nonlinear optics imposes a certain time ordering by singling out one of the fields. This 
response is causal. The scattering description is noncausal since it allows for arbitrary time 
ordering of interactions with the various fields. 
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1 1 
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